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ABSTRACT 
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The thermal condition of cryogenic fluids in the vertical dis- 
charge line during draining of saturated liquid from a self-pressurized 
container was investigated. At the discharge line inlet the fluid be- 
comes two-phase; however, due to the increased static pressure 
while flowing downward in the discharge line the vapor subsequently 
condenses. The per cent of evaporation and the distance required to 
recondense are dependent upon the type of inlet, rate of discharge, 
and initial fluid properties. Equations are derived to predict the 
fluid quality and static pressure in the discharge line. Predictions 
are made for the maximum fluid quality, the distance required for 
vapor condensation, and the vapor variation with elevation. Analyti- 
cal and experimental results are compared for moderate 
pressures and discharge line flowrates. 
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Definition 
Cross Section Area (m 2 ) 

Slope of linearized saturated liquid enthalpy curve 
versus pressure 

Intercept of linearized saturated liquid enthalpy 
curve versus pressure 

Dynamic pressure (N/cm 2 ) 

Discharge line diameter (m) 
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Gravitation constant (9.81 m /sec 2 ) 

Enthalpy (J/G) 

Enthalpy of vaporization (J/G) 

Length (m) 

Static pressure (N/cm ) 

Change in static pressure (N/cm ) 

Volumetric flowrate (cubic meters per minute) 
Heat transferred to fluid (J/G) 
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Specific volume (m 3 /kg) 
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Definition 
Gravimetric quality (%) 

D0 

Parameter defined by Martinelli 
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Parameter defined by Martinelli 
Volumetric quality (%) 

Elevation (m) 

GC 

Parameter defined by Martinelli 
Viscosity (N sec/m 2 ) 

Surface tension (N/cm) 

Velocity (m/ sec) 

SUBSCRIPTS 

Container 
F rictional 
Vapor phase 
Liquid phase 

Arbitrary plane passing through discharge line 
Arbitrary plane passing through discharge line 
Inlet throat 
Two-Phase 

Turbulent liquid -turbulent vapor flow 
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Symbol 

Definition 

vt 

Viscous liquid-turbulent vapor flow 

vv 

Viscous liquid- viscous vapor flow 
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Elevational 


x 


NASA-GEORGE C. MARSHALL SPACE FLIGHT CENTER 


TECHNICAL MEMORANDUM X- 53330 


FLUID QUALITY IN A SELF-PRESSURIZED 
CONTAINER DISCHARGE LINE 

By 

Hugh M. Campbell, Jr. 


SUMMARY 

The behavior and conditions of the fluid in the vertical discharge 
line below a closed tank from which a saturated liquid is being drained 
were investigated. The saturated liquid, upon entering the discharge 
line, becomes two-phase due to the velocity component of pressure 
and subsequently condenses due to increasing liquid head as it travels 
vertically downward in the container discharge line. 

Equations were derived from Bernoulli's Equation, the general 
energy equation, and the equation of state to express the discharge 
line static pressure and fluid quality as a function of container pressure, 
discharge line flowrate, elevation, and fluid properties. Solutions to 
the equations were obtained by a method of finite differences and by an 
analytical method. Frictional and inlet pressure losses were considered 
only in the solutioft obtained from the method of finite differences. 
Experimental results are presented and compared with the analytical 
results . 


INTRODUCTION 

Present day space vehicles require pressurization of propellant 
tanks to insure adequate net positive suction head at the propellant 
pump inlet and to maintain vehicle structural integrity. Adequate net 
positive suction head must be maintained to prevent pump cavitation, 
for pump cavitation causes erratic propellant flow that can result in 


catastrophic failures. Propellant tanks are also pressurized to assist 
the thin tank walls in maintaining structural integrity. A loss of 
propellant tank pressure during powered flight could result in complete 
structural failure. 

Generally, cryogenic propellant tanks are filled at atmospheric 
pressure, and just prior to flight, are pressurized to subcool the 
liquid. By this process, liquid vaporization is prevented, and more 
favorable pump inlet conditions are maintained. Container pressure is 
maintained during draining by the introduction of a pressurant into the 
container from an external source. 

When a superheated vapor below critical pressure coexists in a 
container with the same fluid in a subcooled state, a saturated liquid 
layer will form at the vapor-liquid interface, which is in thermal 
equilibrium with the vapor, and thermal gradient will extend down- 
ward into the liquid. Any disturbance of the vapor-liquid interface 
brings the subcooled liquid into contact with the superheated vapor, 
causing vapor condensation with an accompanying drop in tank pressure. 
This phenomenon has been observed in a container that was oscillated 
to produce sloshing of the liquid and in some propellant piping systems. 
In the Saturn systems an inert non- condensing pre-pressurant is used 
to minimize the effects of the contact between superheated vapor and 
subcooled liquid. 

An alternate pressurization scheme consists of filling the cryo- 
genic container, closing off the container, and allowing the heat that 
is transferred to the cryogen to saturate the cryogen and increase the 
container pressure to the required amount. The container can then be 
drained without the addition of pressurant from an external source. 
Pressurant addition to the ullage space comes from the evaporation of 
the bulk liquid at the vapor-liquid interface. As the container pressure 
drops, sensible heat from the bulk liquid supplies the heat for evapora- 
tion. Thus, during container draining, as heat is continuously removed 
from the bulk liquid, the liquid temperature and tank pressure con- 
tinuously decrease at a slow rate. 

Utilization of self-pressurization in discharging cryogenic pro- 
pellants precludes any pressure drop due to liquid-vapor interface 
disturbance, because the liquid is saturated and the vapor is only 
slightly superheated. No pressurant need be relieved, since the 
highest container pressure obtained would be the pressure required 
initially. 
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The primary disadvantage to the method of self-pressurization 
is the existance of two-phase fluid flow in the container discharge line. 
This reduces the mass capacity of a given line size and complicates 
the accurate control of propellants to the combustion chamber. 

In this report the quality of the cryogenic propellant in the 
vertical discharge line (through which a self-pressurized propellant 
container is being drained) as a function of container pressure, 
discharge line flowrate and location, is investigated including the 
influence of the increased static pressure (liquid head) on a fluid 
element flowing vertically downward in causing the vapor phase to 
condense. 


RELATED INVESTIGATIONS 

In I960 Canty^ analyzed the thermodynamics of a cryogenic 
container, initially full of saturated liquid, when drained adiabatically 
without the addition of a pressurant. Curves were presented depicting 
ullage pressure as a function of the initial pressure and the liquid 
(by weight) removed from the container. A copy of these curves for 
oxygen and nitrogen are shown in FIG 1. In the energy balance, Canty 

assumed a saturated liquid leaving the container. Based on this 
assumption, as the saturated liquid incurs a velocity in leaving the 
container, the velocity or dynamic pressure will increase at the 
expense of static pressure, causing two-phase flow in the discharge line. 
While flowing downward in the discharge line, the fluid static pressure 
increases due to the elevation difference, thus condensing the vapor 
phase if frictional pressure losses are not too great. 

TWO-PHASE FLUID FLOW 

Two-phase fluid flow of a monocomponent fluid is the simultaneous 
flow of a liquid and vapor in a common pipe or tube. Two-phase fluid 
flow exhibits interesting patterns and is frequently classified according 
to pattern and the flow condition of each phase. A survey of the literature 
on two-phase fluid flow in a horizontal tube by Leonhard and McMordie^ 
revealed the following flow patterns. 

1. Bubble Flow: Flow characterized by the formation of 

individual bubbles along the upper surface of the tube. Bubble 
size increases as the fluid flows through the tube. 
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Z. Froth Flow: A foam-like mixture of small bubbles 

intimately mixed with the liquid. Froth flow approaches a 
homogeneous fluid. 

3. Stratified Flow: The flow of the liquid phase along the 

bottom of the tube, while the vapor phase occupies the upper 
portion of the tube. The vapor phase flows at a much higher 
velocity than the liquid at the bottom. The vapor-liquid inter- 
face is a sharp, well defined plane. 

4. Wavy Flow: Similar to stratified flow except that waves 

are formed at the vapor-liquid interface. 

5. Plug Flow: Flow characterized by the movement of large 

plugs of liquid through the tube. Plug flow is the pattern likely 
to exist in the transition from bubble to stratified flow. 

6. Slug Flow: Alternate sections of vapor and liquid flowing 
through the tube. It occurs as the vapor flow increases during 
plug flow. 

7. Annular Flow: The flow of a continuous liquid layer along 

the tube wall. The vapor flows along the central core of the 
tube at a much higher velocity than the liquid. 

8. Mist Flow: The flow in which fine liquid drops are suspended 

by surface tension. There is no significant relative velocity be- 
tween the liquid and vapor phase. 

Gravititational force tends to pull the liquid phase to the bottom 
of a horizontal tube and allows the vapor phase to flow along the top of 
the tube, because the liquid is more dense than the vapor. The gravi- 
tational force then causes wavy, stratified, plug, slug, and bubble 
flow. In a vertical tube the gravitational force does not separate the 
phases but rather causes the buoyant force of the vapor phase to 
oppose downward flow. This causes froth, slug flow, or bubble flow, 
in which the bubbles are evenly distributed within a typical tube cross 
section. A map of two-phase fluid flow pattern regions for a 
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horizontal tube is shown in FIG 2. 



(T H, V g V f W g 

FIG 2 Flow Pattern Regions 


Two-phase fluid flow is also classified into the following regimes; 

1. Liquid phase turbulent, vapor phase turbulent 

2. Liquid phase turbulent, vapor phase viscous 

3. Liquid phase viscous, vapor phase turbulent 

4. Liquid phase viscous, vapor phase viscous 

The criterion for classifying two-phase fluid flow into one of 
these four regimes is the Reynolds number computed as though the 
phase in question occupied the complete tube cross sectional area. 
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TWO-PHASE FRICTIONAL PRESSURE LOSS 

Frictional pressure losses in two-phase fluid flow have been 
studied by many investigators using unique hypotheses and approaches. 
Some of the earlier and more important contributions are discussed 
below. 

McAdams, • Woods and Heromar J analyzed the heat transfer 
coefficient and pressure drop for benzene and lubricating oil flowing 
inside a heated horizontal tube. They assumed equal velocities for 
the liquid and vapor phase with a Fanning expression for wall drag 
and an energy balance to obtain apparent friction factors. 

"4I 

Benjamin and Miller^" analyzed the pressure drop for flashing 
water, using the general energy equation and assuming isentropic 
expansion. The lack of experimental data precludes the use of their 
analysis for anything other than background literature. 

r 5 -> 

Davidson et aF w studied heat transmission and frictional pressure 
drop in boiler tubing using the Fanning frictional pressure loss equa- 
tion and experimentally determining the friction factor in terms of 
Reynolds number. 

r 6 

In 1944, Martinelli et al~ " made a major contribution to two-phase 
fluid flow analysis techniques using the Darcy- Weisbach frictional 
pressure loss equation and assuming equal static pressure losses for 
the liquid and vapor phase. A dimensionless parameter ^g was 
defined so that: 



or 


and 



U) 

(3) 

(4) 
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For turbulent liquid, turbulent vapor flow, the parameter is defined 


L tt 


/ II \ 0111 / v \ 0 555 / 

■ (~j m 

g' ' g' ' g/ 


(5) 


For viscous liquid, turbulent vapor flow, the paramete r Y vt is defined 

4> 2 g versus and Y vt were determined experimentally using air and 
eight assorted liquids, including benzene, water and SAE-40 oil. 

In 1945, Martinelli, Putnam and Lockhart^ presented data 
depicting 4^ in terms of Y vv where Y^ v is defined by 

Yvv = ijt) (*7) (” 7) (7> 

^ g was determined experimentally in terms of Y vv using air and oil 
flowing through a capillary tube. 


In 1948, Martinelli and Nelson® predicted the frictional pressure 
loss for forced circulation of boiling water using the relationship 


tt 


/ v0.143y v0.57L v 

■a &) (4- - ■) 


(8) 


with curves for and versus Ytt from reference 5. The predicted 
and measured pressure losses agreed within j_ 30 per cent in most 
cases . 


CsD 

In 1949, Lockhart and Martinelli u recorrelated the various 
parameters related to frictional pressure losses in two-phase fluid 
flow. In this recorrelation the parameter was defined so that 



Y tt = 

l.n 

x tt 


(9) 

or 

x tt = 

/jM° 1 L 

. 0.5 . . 0.9 

lf_\ /*| \ 

(10) 


\ n g / \ 

v g / \ w g/ 

V s 

plotted versus 

x tt in FIG 

3. 
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Gazely and Bergelen, in a discussion of Lockhart and Martinelli's 
work, stated that they had obtained data from independent tests using 
air and water. This data, when analyzed using the parameters of 
Lockhart and Martinelli, agreed within +20 and -30 per cent of 
Lockhart and Martinelli 1 s curves. 

Throughout Martinelli and associates work a wide spread of data 
exists, indicating the omission of a pertinent parameter. Gazely and 
Bergelen postulate that the pertenent parameter could be velocity. 
Martinelli and associates' data exhibits considerable scatter; however, 
theirs is the best approximation revealed in the literature. Their 
analysis is frequently used by other authors for computing two-phase 
fluid flow frictional pressure loss. 

Johnson and Abou-Sabe ^ studied heat transfer and frictional 
pressure loss for two-phase fluid flow. Their data agreed within +25 
and -50 per cent with the data of Martinelli and associates. 

Harvey and Foust®® mathematically derived a relationship for the 
peak temperature and the vapor head temperature, attained in a long 
vertical tube evaporator, using the equations of state, energy balance, 
momentum balance and equation of continuity. The analysis of 
Martinelli and associates was used to predict frictional pressure loss. 
Experimental verification proved the analysis, "accurately applicable 
to experimental data for water below atmospheric pressure". 

Rogers^? analyzed the flow of flashing liquid hydrogen based on 
the data and analysis of Harvey and Foust, and Lockhart and Martinelli 
No experimental data were presented; however, the following equations 
were derived using the data of Lockhart and Martinelli. 

<t>J= 13.054X 1 t J 6 + 2.5996 for 0SX<0.9 (11) 

4>*= 15. 843X°' t 659 for 0.9^ X<26.0 (12) 

Hatch and Jacobs'- 1 ® studied frictional pressure loss in two-phase 
fluid flow using trichloromonofluoromethane and hydrogen. Their data 
was 10 to 30 per cent lower than that of Martinelli and associates when 
4* i versus X was plotted. Hatch and Jacobs attribute part of this incon 
sistency to the probable existence of thermodynamic metastability. 
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CRYOGEN PROPERTIES 


The physical and thermodynamic properties of many of the cryo- 
genic fluids and the properties of many structural materials at cryo- 
genic temperatures have been studied extensively by the National 
Bureau of Standards, Cryogenic Engineering Laboratories. 

r i41 

Johnson published curves showing density, expansivity, thermal 
conductivity, specific heat, enthalpy, heat of vaporization, heat of 
fusion, heat of sublimation, phase equilibria diagrams, dielectric 
constants, ads orptivity, surface tension and viscosity of most of the 
cryogenic fluids. The thermal expansion, thermal conductivity, 
specific heat and enthalpy of many engineering materials were also 
shown. Viscosity, from this reference, for saturated liquid nitrogen 
and oxygen is plotted versus pressure in FIG 4. 

jj 1 

Stewart, Hust and McCarty J tabulated the enthalpy, density and 
entropy for liquid and gaseous oxygen. The enthalpy versus pressure 
for saturated liquid oxygen from the reference is plotted in FIG 5, 
enthalpy of vaporization versus pressure in FIG 6, specific volume 
for saturated liquid versus pressure in FIG 7, and specific volume 
for saturated vapor versus pressure is shown in FIG 8. 

Strob ridge tabulated the enthalpy, internal energy, entropy and 
specific volume of liquid and gaseous nitrogen. From this reference, 
enthalpy for saturated liquid nitrogen is plotted versus pressure in 
FIG 5, enthalpy of vaporization versus pressure in FIG 6, specific 
volume for saturated liquid versus pressure in FIG 7 and specific 
volume for saturated nitrogen vapor is shown plotted versus pressure 
in FIG 8. 


ANALYSIS 


ANALYTICAL SOLUTION 

Fluid Conditions in Pipe Inlet 

Consider a self pressurized cryogenic container from which a 
saturated liquid is being drained through a vertical discharge line. 

As the fluid flows past the container boundary (see FIG 9) it is assumed 
to be in the liquid state with negligible dynamic pressure. Additional 
assumptions are (a) thermodynamic equilibrium, (b) homogeneous fluid 
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FIG 5 Enthalpy versus Pressure for Saturated Liquid Nitrogen and Oxygen 
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FIG 6 Enthalpy of Vaporization versus Pressure for Nitrogen and Oxyg 
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FIG 7 Specific Volume versus Pressure For Saturated Liquid Nitrogen and Oxygen 
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properties (lumped parameters), (c) no heat transfer to or from the 
fluid in the inlet, (d) negligible elevation effects and (e) no frictional 
or inlet pressure losses. The saturated liquid enthalpy is assumed 
to be a linear pressure function over short pressure ranges, and the 



FIG 9 Container, Pipe and Inlet Boundaries 


liquid and vapor specific volume and enthalpy of vaporization are 
assumed to be constant over short pressure ranges. The coordinate 
system is so chosen that z is positive downward with the origin at the 
inlet throat. 

The differential form of Bernoulli's equation is 


-4 -4 

, 1 0 10 gdz 

ap + - — - + s = U 

v v 


(13) 


When elevation effects are excluded, equation 13 becomes 


dp + iolIm = o 


(14) 


17 


The equation of continuity is 
^ = (w/A) v 


( 15 ) 


In differential form equation 15 becomes 

d4* = (w/A)dv+ vd(w/A) (16) 

Substituting equation 15 and 16 into equation 14 and rearranging, 

10 4 dp + (^-fdv+ v(-g-)d(-^-) = 0 (17) 

Equation 17 is subject to the boundary conditions: 

P = P c 

x = 0 

(w/A c ) = 0 


The specific volume of a two-phase fluid in terms of saturated liquid 
and vapor specific volume and gravimetric quality is 


v = ( 1 — x) v. + 


xv 


g 


(18) 


In differential form, considering constant liquid and vapor specific 
volume, equation 18 becomes 

dv - (v - v . ) dx 

g i 

Substituting equation 18 and 19 into equation 17 gives 

10 dp + (w/A) (Vg - v^)dx + J(Vg - v^)x + v^j (w/ A)d(w/ A) = 0 
In differential form the general energy equati< 


ion is 


-3 -3 

dh + 10 ^d^ + 10 gdz + dg - dW = 0 


(19) 


( 20 ) 


( 21 ) 


Using the assumptions stated at the first of this chapter, equation 21 
b acornes 


dh = 0 


( 22 ) 
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The enthalpy of a two-phase mixture in terms of saturated liquid 
enthalpy, enthalpy of vaporization, and gravimetric quality is given by 


h = h # + xAh 


(23) 


But 


So 


£ = ap + b 


h = ap + b + xAh 

In differential form equation 25 is 

dh = adp + Ahdx = 0 

To determine x^. equation 26 is solved for dp and this value is 
substituted into equation 20. Equation 20 then becomes, after re- 
arranging. 


(24) 

(25) 

( 26 ) 




dx + 


(v g - v t ) x + v ( 


( ir )d( a~ ) = 0 (27) 


Solving equation 27 by separation of variable^ the result is 

P , ^ , A l '/ 2 

[Vg - v ,) x+ v f 


2 4 

(-~-) (v„ - v.) - U Constant (28) 

A&* a J 


After incorporation of the boundary conditions, equation 28 becomes 


*t 


= 11 

V V i 


10 4 Ah 


(w/A) (vg-Vj) -lO Ah 


1/2 


v f 


v g -v £ 


(29) 


To determine p^, equation 26 is integrated subject to the boundary 
conditions 

x = 0 


P = P, 


The result is 


Xj. Ah 


P t = P c - a 


(30) 
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Substituting equation 29 into equation 30 gives 


Pt = Pc + 


Ahvj 



v f )a 


Ahv > 

(Vg - v ( )a 


-10 4 Ah/a 

(w/A) 2 (v g - v £ )-10 4 Ah/a 


(31) 


The fluid conditions in the pipe inlet throat are described by equations 
29 and 31. These equations will be used with equations developed in 
the next section to determine the fluid conditions in the cylindrical pipe. 


Fluid Conditions in Vertical Pipe Below Inlet 


The fluid will cross the exit boundary of the inlet with a quality 
X|- and pressure p^ determined using equations 29 and 31. Assuming 
a cylindrical pipe and the same assumptions listed in the previous 
section, except the velocity effects are neglected and elevation effects 
are included in Bernoulli's equation, equation 13 becomes: 


dp + 


- 4 

10 g dz 
v 


= 0 


(32) 


Equation 32 is subject to the boundary conditions 


z = 0 
x = x+ 


The specific volume of a two-phase mixture in terms of saturated 
liquid and vapor specific volume and gravimetric quality is 


v = (Vg - Vj ) X + v £ 


(33) 


Substituting equation 33 into equation 32 and rearranging, 

10 4 (Vg - Vj)x + v,] dp + gdz = 0 (34) 

For the case in question, the general energy equation reduces to 

dh = 0 (35) 
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Using a linearized relationship for the saturated liquid enthalpy over 
short pressure ranges, substituting this into the relationship for the 
enthalpy of a two-phase mixture, differentiating, and solving for dp 
yields: 


dp = - 


Ahdx 

a 


(36) 


Substituting equation 36 into equation 34 and rearranging 

— — |^( v g " x + dx ~ dz (37) 

Integrating equation 37 and incorporating the boundary conditions 

2 

z = 1 ^ir L { K • v « )Xt + v <] ' [ (v g • v ( )x + 'i] } (38) 

To determine the distance required to condense the vapor, set x = 0 
in equation 38, which yields 


10 4 Ah i\, N , 

= -lir i r« ■ v * )Xt + v < 




(39) 


Solving equation 38 for x, the quality variation along the discharge 
line is given by 


x = 


v g - v f 


(v g - V + y t 


(2 x 10 4 )gaz 


Ah 


0.5 


v t 


V g- V t 


(40) 


Integrate equation 36 to determine the variation of quality with pressure; 
then, using the boundary conditions, 

x* Ah Ahx 


P = Pt 


(41) 


Substituting equation 40 into equation 41 

x^ Ah Ahvj 


P = Pt 


a(v. 


( v g - v | )x t + v ( 


j (2 x 10 )gaz j | + 5 _ 


Ahv* 


(42) 


a(v g - v f ) 
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By using equations 29 and 31 with equations 40 and 42, the fluid 
conditions can be determined at any discharge line elevation in terms 
of elevation, container pressure, and weight flowrate per unit cross 
sectional area. 

FINITE DIFFERENCE SOLUTION 
Fluid Conditions in Pipe Inlet 

Consider a container from which a saturated liquid is being 
drained through a vertical discharge line. As the fluid flows past the 
container boundary (FIG 10), it is assumed to be in the liquid state 
with negligible dynamic head (pressure). Additional assumptions are 
(a) thermodynamic equilibrium, (b) homogeneous fluid properties and 
(c) no heat transfer to or from the fluid in the pipe. The coordinate 



FIG 1.0 Container, Pipe and Inlet Boundaries 

system is so chosen that the positive direction is vertically downward 
and the origin is at the inlet throat. The pressure of the fluid, as it 
crosses the container boundary, can be related to the pressure at the 
inlet throat by 


Pc + A Pz = Pt + ^ + A Pf 


(43) 


22 



No literature information has been available to predict the frictional 
and inlet pressure losses between the container and throat, however, 
in Reference 17, a value of 4 per cent of the dynamic head was recom- 
mended for a well shaped bell inlet. From this recommendation, 
frictional and inlet pressure losses were assumed to be 10 per cent 
of the dynamic head or. 


-4 2 

1.1x10 4*t 

Dt + Apr = 

v 


(44) 


Combining the equation of continuity and the relationship for specific 
volume of a two-phase fluid in terms of saturated liquid and vapor 
specific volume and gravimetric quality gives 


D 


+ A P f = (1.1 x 1 0 ) (w/A) 


< 1 - x t > v n + h Vj 


(45) 


The pressure gain due to elevation is given by 

z„ 


Ap 7 = g x 1 




c 

dz/ v 


(46) 


For small differences in specific volume and elevation, the pressure 
gain due to elevation can be approximated by: 

(g x 10 " 4 ) [z t - zj 


A P ; 


0. 5 


v + v 
c t 


(47) 


Substituting the relationship for specific volume of a two-phase 
mixture in terms of saturated liquid and vapor specific volume and 
gravimetric quality, noting that the quality crossing the container 
boundary is zero, gives 


-4 


A Pz = 


( 2 gx 10 ) [z t - z c ] 


(48) 


v tc + (1 - X t )v /t + x t v gt 
Substituting equations 45 and 48 into equation 43 and rearranging gives 


(2g x 10 ) f z . - z 1 _ 4 , j 

p t = P c + — r: ; - (1.1 X 10 ) (w/A) 


v fc +(l-x t )v it +x t v 


gt 


(1 - x t )v| t + x t v 


gt 


(49) 


The general energy equation, as applied to the pipe inlet, is 



Using the assumptions listed, with the relationship for the enthalpy 
of a two-phase mixture in terms of saturated liquid enthalpy, enthalpy 
of vaporization and gravimetric quality, and noting that the quality of 
the fluid crossing the container boundary is zero, equation 50 reduces 
to 



(51) 


(52) 


Equations 49 and 52 were solved on a digital computer using true 
fluid properties. Specifying the container pressure and discharge line 
weight flowrate per unit cross sectional area, and assuming a throat 
quality x t , equation 49 is solved for throat pressure. This throat 
pressure is then used in equation 52 to determine a corrected qualtiy. 

The corrected throat quality is used in equation 49 to determine a 
corrected throat pressure. This process is repeated until two con- 
secutive pressures and qualities agree within 0. 1 per cent of their 
absolute value. Throat quality is plotted versus container pressure 
for several values of w/A for nitrogen in FIG 11 and for oxygen in FIG 12. 


Upon examination of equation 52, noting that the enthalpy of a 
saturated liquid is a function of pressure, it can be seen that the 
throat quality will be zero, but the fluid will be saturated when the 
static pressure in the throat and container are equal. Substituting 
these conditions into equation 49, and noting that since the throat and 
container pressures are equal, their specific volumes will be equal, 
and the following relationship is obtained. 


w/A 


,[ 8.918 (z t - z c )] °- 5 
v *c 


(53) 


In the experimental facility used to verify this analysis, the term 
z t “ z c was 0.254 meters. Using the value in equation 53, the 
maximum value of w/A, for which the fluid will remain single phase, 
is plotted versus container pressure in FIG 13. Any value of w/A 
greater than the plotted value will result in two-phase flow, and any 
value below the plotted value will cause subcooled fluid in the pipe 
inlet throat. 
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Pressure (N/cm 2 ) 

FIG 11 Nitrogen Inlet Quality versus Container Pressure (Finite Difference Solution) 
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FIG 13 Maximum Single Phase w/A in Pipe Inlet versus Container Pressure For Nitrogen and Oxygen 



Fluid Conditions in Cylindrical Pipe Below Inlet 


The fluid crosses the exit boundary of the inlet with a quality and 
pressure that can be determined by using equations 49 and 52. The 
coordinate system is so chosen that the positive direction is vertically 
downward and the origin is the inlet throat. Using the same assump- 
tions as in the previous section, the fluid pressure at the arbitrary 

elevations, z and z , where z > z can be related bv 
m n m n 3 


Pm ^m A Pf Pn ^n Ap. 


(54) 


where 


D = 


10"V 


(55) 


Using the equation of continuity and the relationship 
v = (1 - x)vj + xv g 
equation 55 becomes 


D m = 10 ( W / A ) 2 " X ~J v #_ + v . 


^ ) Vg T X V 

m ' m m gm 


-4 


D n = 1° ( w/A ) [f 1 - x n) v fn + x n 

The pressure gain due to elevation is 


g n 


(56) 


(57) 


Ap = 1 


_4 r m 
0 g I dz/v 


n 


(58) 


For small elevations and specific volume changes, the pressure gain 
due to elevation can be approximated by 


Ap = 


2g xl ° (z m - Zn ) 


V + V 
m T n 


(59) 


Substituting equation 56 into equation 59 


A P.. = 


( z m z n) ^0 


-4 


(1 - X )v. + x V + (1 - X )v, + X V 

m l m m gm n' in n gn 


(60) 
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To determine the frictional pressure loss, the analysis of Martinelli 
and associates 9 > 9 is used. All discharge line flow is assumed to 
be turbulent liquid, turbulent vapor, and only the equations for this 
case are selected. From reference 6 the relationship for single 
phase and two-phase frictional pressure loss is 


dp 

dz t P 


dp 

dz 



( 61 ) 


is determined from the following equations: 

-1. 36 

= 13. 054 X + 2.5996 for 0iX<0.9 


where 


r t 


4>J = 15.843 X 


-0. 659 


dj = 4.957 X 
4> 2 t =l 


-0. 303 


X = 


\°- s , . . 0.1 . , . 0.9 

1 1 




for 0 - X < 26 . 0 

for 26.0 <X<196.0 
for X 2 196 . 0 


>(62) 


(63) 


Since v^, v , p^and p„ are functions of pressure for a saturated 
fluid, the function (vf/' Vg)°‘ 5 ( Pj/ P g) 0- ' is also a function of pressure, 
and is plotted versus pressure in FIG 14. 


The Darcy- Weisbach relationship for frictional pressure loss is 

. , 
dp 

dz 


4f + (10 ) 

l v l d 

The Blasius relation shin for the friction factor for turbulent flow is 


(64) 


f = 


0. 046 


R = 


R 

wd 


0.25 


Ap 


(65) 

( 66 ) 


t 


Substituting the equation of continuity, equation 65 and 66 into equation 
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s - 0 (*a/®a) 


64, then sutstituting this into equation 61 gives 


dp 

dz 


0. 184 (w/A^Vj (10~ 4 ) <j>£ 


tp (wd/ A(i| ) 0,25 d 


( 67 ) 


The frictional pressure loss between elevations z and z is 

n m 


^ -f-t 


dz 


( 68 ) 


n 


tp 


For small changes in frictional pressure loss and elevation, the 
frictional pressure loss can be approximated by 


Apf -0.5 (z m - z n ) 


dp 

, dp 

! dz 

+ ~Kz~\ 

tpm QZ > tpn 


(69) 


Substituting equation 67 into equation 69 and rearranging. 


_ 9.2(10~ 6 ) (z m - z n ) (v /m fy 2 m + v| n 4>| n ) (w/A)' 

n 95 

(wd/Ap| ) ‘ d 


(70) 


Substituting equations 57, 60, and 70 into equation 54 and rearranging^ 

p = p + (w/A) 2 (1-x )v- + x v - (1 - x )v # - x v 

^m ^ n [j m In n gn ' m l m m gm 


-4 


2g( z m " z n) ( 10 ) 


(1 - X )Vj + X V + (1 - X )v- + X 

v rrv t m m gm n In 


n V gnJ 


(71) 


9. 2(1 0“ 6 ) (w/Af(z m - z n ) (v /m 4>ij m + v /n 4> 2 , n ) 


. 0.25 

(wd/ Ap ^ ) d 

The general energy equation relating the energy flowing past, the 
container boundary to the energy flowing past the elevation plane z m , 
after incorporation of the assumptions listed earlier, is 


^c 


Substituting equation 51 into equation 72 and solving for x m gives 

k/c " k-im 


x 


m 


Ah. 


(72) 


(73) 


m 
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Using solutions obtained by equations 49 and 52 as a starting 
point, the fluid conditions were computed at discrete discharge line 
elevations in the same manner that equations 49 and 52 were solved. 
For all combinations computed of container pressure and discharge 
line flowrate per unit of cross sectional area, the quality was shown 
to decrease as the fluid moved down the discharge line. The distance 
required to condense the vapor is plotted versus container pressure 
for several values of (w/A) for nitrogen in FIG 15 and for oxygen in 
FIG 16. The variation of quality with elevation is shown for two con- 
tainer pressures and several values of (w/A) for Nitrogen in FIG 17 
and for oxygen in FIG 18. 

EXPERIMENTAL APPARATUS AND PROCEDURE 
EXPERIMENTAL APPARATUS 
General Description 

The experimental apparatus consisted of a 1.702 meter (68 inches) 
inside diameter, spherical aluminum tank atop a 12.2 meter (40 foot) 
tower. A 0. 1524 meter (6 inches) inside diameter discharge line and 
a 0. 1015 meter (4 inches) inside diameter recirculation line were 
suspended vertically underneath the tank and connected near the 
bottom of the discharge line. The discharge line and inlet were 
insulated with polyurethane foam to minimize heat transfer. The 
recirculation line was not insulated, but heat was allowed to flow into 
the fluid so that a density differential would exist between the fluid 
inside the recirculation line and the fluid inside the discharge line to 
promote convective circulation and prevent "geysering". The 
experimental apparatus is depicted schematically in FIG 19, and a 
picture of the apparatus is shown in FIG 20. A recirculation valve was 
used to prevent flow in the recirculation line during draining of the 
tank. A vent valve was provided atop the tank to relieve tank pressure, 
and a valve at the bottom of the discharge line was used to drain the 
container. Discharge line flow was controlled with orifices located 
immediately upstream of the drain valve. A description of the com- 
ponents used in this apparatus is shown in Table I. 

Pressure was measured at two locations in the container, and 
temperature was measured at five locations. Liquid depth was 
determined by measuring the difference in static pressure at the top 
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Pressure (N/cm 2 ) 


16 Distance Required to Condense Oxygen Vapor versus Container Pressure 
(Finite Difference Solution) 








Elevation (m) 


17 Quality versus Elevation for Nitrogen (Finite Difference Solution) 





















P -X Denotes Pressure 
Measurement 
Q-X Denotes Quality 
Measurement 

T-X Denotes Temperature 
Measurement 


Vent Valve 


Tank I 


T-2C 


T-4 <T 
T-5<C 


L» lL i nn 


1 nn 

AU Inlet Throat (Z=0) 


Recirculation Line 


Recirculation Valve 


Discharge Line 


Cryogen 

Supply 


Fill Valve 


P-10 


O O | Flowmeter 


Orifice 



Drain Valve 







Table 
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Table II 

Measurement Locations 

Measurement 

Measured 

Location 

Number 

Property 


F-l 

Discharge Line 
Flowrate 

In discharge line 11.66 m below inlet throat 

P-1 

Static Pressure 

On top of tank measuring ullage pressure 

P-2 

Static Pressure 

On tank bottom measuring container pressure 

P-3 

Static Pressure 

On discharge line 0. 076 m below inlet throat 

P-4 

Difference in 
Static Pressure 

On discharge line between 0.076 and 0.908 m 
below throat 

P-5 

Static Pressure 

On discharge line 2.228 m below inlet throat 

P-6 

Difference in 
Static Pressure 

On discharge line between 2.228 and 3.060 m 
below throat 

P-7 

Static Pressure 

On discharge line 4.383 m below inlet throat 

00 

1 

Difference in 
Static Pressure 

On discharge line between 4.383 and 5.212 m 
below throat 

P-9 

Static Pressure 

On discharge line 7.861 m below inlet throat 

P-10 

Static pressure 

On discharge line 10.71 m below inlet throat 

P-11 

Difference in 
Static Pressure 

Between top and bottom of tank 

Q-l 

Fluid Quality 

In discharge line 0.457 m below inlet throat 

Q-2 

Fluid Quality 

In discharge line 2.61 m below inlet throat 

Q-3 

Fluid Quality 

In discharge line 4. 76 m below inlet throat 

T-l 

Temperature 

In tank on vertical axis 0.762 m above tank center 

T-2 

Temperature 

In tank on vertical axis 0.381 m above geometric 
center 

T-3 

Temperature 

In tank at geometric center 

i 

H 

Temperature 

In tank on vertical axis 0.381 m below geometric 
center 

T-5 

Temperature 

In tank on vertical axis 0.762 m below geometric 
center 

T-6 

Temperature 

In discharge line 0.076 m below inlet throat 

T-7 

Temperature 

In discharge line 0.908 m below inlet throat 

00 

i 

H 

Temperature 

In discharge line 3.060 m below inlet throat 

T-9 

Temperature 

In discharge line 5.212 m below inlet throat 

T- 10 

Temperature 

In discharge line 7.861 m below inlet throat 

T-ll 

Temperature 

In discharge line 10.711 m below inlet throat 
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and bottom of the container. The temperatures inside the tank were 
measured with triple junction copper constantan thermopiles located 
symmetrically along the vertical axis of the tank, a distance of. 0.381 
meters (15 inches) apart. The thermopiles were mounted on a 
stainless steel bracket suspended from the top of the tank, and the 
wires were led to a reference temperature junction located outside 
the tank through an adapting flange between the tank and vent valve. 
Pressure was measured at eight locations in the discharge line. 
Temperatures were measured along the center line of the discharge 
line at eight elevations, and fluid quality was measured at three 
discharge line elevations. The location of all discharge line measure- 
ments is shown in Table II as a vertical distance below the inlet 
throat. All measurements were recorded on a 36- channel oscillograph. 

A discussion of the measurement equipment used in recording data from 
this experimental apparatus follows. 

Flow Measurement 

The discharge line flowrate was measured with a turbine type 
flowmeter. The location of the flowmeter and experimental conditions 
were so chosen that single-phase flow would exist at the flowmeter. 

Pressure Measurements 

Two types of pressure transducers were used for pressure 
measurements; one type (used for measuring P-4, P-6, P-8 and P-11) 
used a twisted Bourdon tube and the other (used for all remaining 
pressure measurements) was a diaphram type. 

Quality Measurements 

The liquid vapor quality was determined by measuring the difference 
in dielectric constant of the liquid and vapor. The transducers were 
pipe test sections with an internal capacitive matrix and temperature 
sensor. The matrix measured the capacitance of the fluid passing 
through the pipe test section. This fluid capacitance was then fed 
into an analog computer where it was compared to the capacitance of 
the liquid and vapor. The capacitance of the liquid and vapor were 

generated as a function of temperature using the output of the tempera- 
ture sensor in the pipe test section. The volumetric liquid vapor quality 
appeared as a direct current output that was recorded on the oscillograph. 
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A block diagram of a typical quality measurement taken from 
reference 18 is shown in FIG 21. 

Temperature Measurements 

All temperatures were measured with triple junction copper 
constantan thermopiles. Liquid nitrogen at atmospheric pressure was 
used as the reference junction. 

EXPERIMENTAL PROCEDURE 

Before conducting an experiment, the measurement systems 
were calibrate^ and the correct orifice was installed in the discharge 
line. The tank atop the tower was then filled from the cryogen 
storage supply, and the system was allowed to remain at rest for 30 
minutes for thermal stabilization. During fill and thermal stabilization, 
the recirculation valve (FIG 19) was open to allow convective circulation, 
and the vent valve* was open to prevent pressure buildup in the tank. 

After termal stabilization, the tank was replenished with liquid 
nitrogen and the vent valve was closed. The container was then 
allowed to pressurize by the flow of heat from the ambient surroundings 
to the fluid through the tank wall and recirculation line. The thermopile 
in the bottom of the tank (T-5, FIG 19) was monitered to insure that 
the liquid temperature was near saturation. When the desired container 
pressure was reached, the liquid was assumed to be near saturation 
temperature; the recirculation valve was closed, recording equipment 
was started, and the drain valve was opened. When the liquid was 
depleted, the drain valve was closed; the recirculation and vent valves 
were opened; recording equipment was stopped; and the residual liquid 
was drained into the cryogen storage supply. The vent valve was open 
until the apparatus returned to approximately ambient temperature, 
then was closed. 


COMPARISON OF EXPERIMENTAL AND ANALYTICAL RESULTS 
EXPERIMENTAL CONDITIONS AND RESULTS 

Three experiments were conducted to verify the analysis. Experi- 
mental parameters were selected so that a range of container pressures 
and discharge line flowrates were considered. The experimental con- 
ditions are shown in Table III, and results are shown in FIGS 22 through 
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FIG 21 Typical Quality Measurement Diagram 











35. The only results presented are for that portion of the discharge 



TABLE III EXPERIMENTAL 

CONDITIONS 


Experiment 

Number 

Orifice 

Diameter 

(cm) 

Initial 
Container 
Pressure 
(n/cm )- 

Barometric 
Pressure 
(mm Hg) 

Fluid 

303001 

7.567 

51.6 

745.5 

Nitrogen 

303002 

8. 158 

36. 3 

760.5 

Nitrogen 

303003 

8. 158 

64.7 

768.6 

Nitrogen 


line in which two-phase flow exists. Conditions in the portion of the 
discharge line in which single-phase flow exists can be accurately 
determined using Bernoulli's equation, assuming constant velocity 
and the Darcy- Weisbach relationship for frictional pressure loss in 
a single-phase fluid. 

DATA REDUCTION 

After reducing data from an experiment, the discharge line flow- 
rate was converted to a weight flowrate per unit cross sectional area 
using the equation 

(w/A) = 737.7 Q (74) 

The value given by equation 74 and the container pressure P-2 (FIG 19) 
were then used in equations 49, 52, 71 and 73 to determine the fluid 
conditions in the discharge line. Pressure from the solution of 
equations 49, 52, 71 and 73 at an elevation of 0.076 meters (3 inches) 
was plotted versus time. This is shown as the dashed line in FIG 22, 
27, and 32. The measured pressure at the same discharge line 
elevation, P-3 (FIG 19), was then plotted and is shown as the circles 
in the same figures. The difference in static pressure between the 
discharge line elevations, 0.076 meters (3 inches) and 0.908 meters 
(35.75 inches), was determined from the solution of equations 49, 52, 
71 and 73 and plotted versus time in FIG 25, 30 and 35. The measured 
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FIG 26 Container Temperature versus Time, Experiment No. 303001 
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values of P-4 (FIG 19) were then plotted for comparison of the 
measured and predicted pressure difference across the capacitive 
matrix. No provisions were made for predicting the pressure loss 
caused by the capacitive matrix. 


Volumetric quality from the quality measurement system was 
converted to a gravimetric quality using the equation 


x 


Jy/vg)_ 


(i - y)l 


v r 


+ (y/v g ) 


( 75 ) 


Equation 75 is derived using a dimensional approach in reference 19. 
Gravimetric quality obtained from the solution of equations 49, 52, 71 
and 73 at a discharge line elevation of 0.457 meters (18 inches) is 
plotted atld compared with those measured at Q-l (FIG 19) in FIG 23, 
28 and 33. 


The discharge line weight flowrate per unit cross sectional area 
is plotted versus time in FIG 24, 29 and 34. 

For additional information on self-pressurized containers, the 
temperatures inside the container are plotted for two of the three 
experiments in FIG 26 and 31. 


CONCLUSIONS 

From comparison of the predicted and experimental results, it 
is concluded that the analytical equations will predict the fluid con- 
ditions in the discharge line when draining a self-pressurized container. 

The quality measurements show that a slight amount of vapor was 
discharged into the discharge line during experiment No. 303001 
between 30 and 40 seconds and during experiment No. 303002 between 
10 and 15 seconds. A subcooled liquid was drained for the first ten 
seconds of experiments No. 3030U1 and 303002. This is evidenced by 
experimental qualities being less than predicted and by the large 
negative slope of the container pressure curve during this interval. 

The frictional pressure drop across the capacitive matrix was expected 
to be large enough to cause slightly higher than predicted qualities. 
These expectations were proven true by the experimental qualities and 
by the difference between the predicted and measured pressure differ- 
ence across the quality measurement matrix. This excessive frictional 
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pressure loss greatly amplified the effect of the vapor discharge into 
the discharge line during experiment No. 303001 and 303002. The 
frictional pressure drop across the capacitive matrix during experi- 
ment No. 303003 was so great that the pressure below the matrix 
(FIG 35) was less than the pressure above the matrix. The measured 
and predicted discharge line static pressure (FIG 22, 27 and 32) 
agreed within the limits of accuracy of the measurement system. The 
container temperatures did not differ significantly from the saturation 
temperature. The experimental phase of this program proved the 
analysis accurate for the case considered, but points out that any 
obstruction in the discharge line will cause a significant deviation 
from the predicted fluid conditions. To obtain extreme accuracy, the 
effects of obstructions must be considered. The assumption of no 
vapor being discharged into the discharge line was valid in almost all 
instances. When vapor was expelled into the discharge line, the 
amount was small and was neglected with a slight loss of accuracy. 

This analytical technique is applicable to many other problems 
involving two-phase fluid flow, if the quality can be predicted at the 
inlet throat. An example of a problem solvable by this technique is a 
liquid discharge line carrying a two-phase mixture because of a 
vortex in the reservoir above the discharge line. In this instance the 
inlet quality can be predicted knowing the liquid drain rate and esti- 
mating the vapor drain rate from the size of the vortex. 
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